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. ABSTRACT 

o ■ 

^ Aims. A detailed study of the spectrum and variability of the source HESS J1745— 290 in the Galactic Center (GC) region using 
i-H new data from the H.E.S.S. array of Cherenkov telescopes is presented. Flaring activity and quasi periodic oscillations (QPO) of 
>^ HESS J1745-290 are investigated. 

5_j Methods. The image analysis is performed with a combination of a semi-analytical shower model and the statistical moment-based 
Hillas technique. The spectrum and lightcurves of HESS J1745— 290 are derived with a likelihood method based on a spectral 
shape hypothesis. Rayleigh tests and Fourier analysis of the H.E.S.S. GC signal are used to study the periodicity of the source. 
Results. With a three-fold increase in statistics compared to previous work, a deviation from a simple power law spectrum is 
detected for the first time. The measured energy spectrum over the three years 2004, 2005 and 2006 of data taking is compatible 
with both a power law spectrum with an exponential cut-off and a broken power law spectrum. The curvature of the energy 
spectrum is likely to be intrinsic to the photon source, as opposed to effects of interstellar absorption. The power law spectrum 
with an exponential cut-off is characterized by a photon index of 2.10 ± 0.04stat ± O.lOsyst and a cut-off energy at 15.7 ± 3.4stat 
± 2.5syst TeV. The broken power law spectrum exhibits spectral indices of 2.02 ± O.OSstat ± O.lOsyst and 2.63 ± 0.14stat ± O.lOsyst 
with a break energy at 2.57 ± 0.19stat ± 0.44syst TeV. No significant flux variation is found. Increases in the 7-ray flux of HESS 
J1745— 290 by at least a factor of two would be required for a 3a detection of a flare with time scales of an hour. Investigation of 
possible QPO activity at periods claimed to be detected in X-rays does not show any periodicities in the H.E.S.S. signal. 

Key words. Galaxy: center - Gamma rays: observations 



Send offprint requests to: Matthieu Vivier e-mail: 
matthieu.vivierOcea.fr, Christopher van Eldik e-mail: 
Christopher . van . Eldik@mpi-hd . mpg . de 

* supported by CAPES Foundation, Ministry of Education of 
Brazil 



1. Introduction 

The discovery of very high energy (VHE) 7-rays from the 
Galactic Cent er (G C) has been reported b y CAN GAROQ 
(Tsuchiya et al l20D4|l . VERITAS ( Kosack et al. [2004)) . H.E.S.S. 
(HESS J1745-290, Aharonian et al. l2004b|l and MAGIC (Albert 
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et al. l2UU6p . Possible associations of the source with the Sgr 
A East supernova remnant (Crocker et al. I2005P and more re- 
cently with the newly detected plerion G359.95— 0.04 (Wang et 
al. I2006P have been widely discussed in the literature. However, 
with the reduced systematic pointing error obtained using 
H.E.S.S. data up to 2006 (van Eldik et al. [2007)) . Sgr A East 
is now ruled out as being associated with the VHE emission of 
HESS J1745-290. The interpretation of the GC TeV signal as 
annihilation products of dark matter (DM) particles has been 
discussed in Aharonian et al. (2006c). It is unlikely that the 
bulk of the signal comes from DM annihilations. One possibility 
is that the supermassive black hole Sgr A* located at the center 
of the Milky Way is responsible for the VHE emission of the de- 
tected HESS J1745— 290 source. In this cut-off in the high 
energy part of the spectrum might be expected (Ballantyne et al. 
I2007i Aharonian and Neronov I2005ap . Time variability as seen 
in X-rays and IR might also appear, along with flares, and QPO 
frequencies such as those found in Aschenbach et al. (|2004|l . 
QPO periods of 100 s, 219 s, 700 s, 1150 s and 2250 s have 
been observed simultaneously in two different datasets collected 
with the Chandra (Baganoff et al. [2002]) and the XMM-Newton 
(Porquet et al. I2003P observatories in 2000 and 2002, respec- 
tively. The power density spectra of the 2003 infrared flare also 
shows three distinct peaks at ~ 214 s, 733 s and 1026 s (Genzel 
et al. l2003ll . fully consistent with three of the five X-ray detected 
periods. Recently however, the validity of the detection of these 
periods has been disputed by infrared observations on the Keck 
II telescope (see Meyer et al. (|2008p for references). 
The search for a curvature in the TeV energy spectrum and 
time variability in the TeV signal is strongly motivated by GC 
wideband emission models and multi-wavelength data, respec- 
tively. In this paper, updated results on the energy spectrum 
and variability of HESS J1745— 290 are presented, based on a 
dataset collected in 2004, 2005 and 2006. Preliminary results on 
the source position and morphology were published in van Eldik 
et al. IP007|) . 



2. H.E.S.S. observations and analysis 

2.1. The H.E.S.S. instrument 

The H.E.S.S. (High Energy Stereoscopic System) instrument 
(Hofmann et al. [2003) consists of four Imaging Atmospheric 
Cherenkov Telescopes (lACTs) located in the Khomas Highland 
of Namibia at an altitude of 1800 m above sea level. H.E.S.S. 
is dedicated to very high energy 7-ray astronomy (defined as 
E > 100 GeV), beyond the energy range accessible to satellite- 
based detectors. The lACTs have a large mirror area of 107 , 
refiecting the Cherenkov light emitted by 7-induced air show- 
ers onto a camera of 960 photomultiplier tubes (PMTs). Each 
PMT covers a field of view of 0.16°. The total field of view of 
the camera is 5°in diameter. The telescopes are positioned at the 
corners of a square, of side 120 m, which allows for an accurate 
reconstruction of the direction and energy of the 7-rays using 
the stereoscopic technique. The energy threshold of H.E.S.S. is 
approximately 100 GeV at zenith. More details on the H.E.S.S. 
instrument can be found in Aharonian et al. (I2006dp . 

2.2. Dataset and analysis 

The previously published H.E.S.S. observations on the GC 
(Aharonian et al. 12006c [) was on a 48.7 hours (live time) data 
sample collected in 2004. In this paper new results on the spec- 
tral analysis of the GC are derived, using data from subsequent 
observation campaigns carried out in 2005 and 2006, with zenith 
angles ranging up to 70°. The runs taken at high zenith an- 
gles are sensitive to higher 7-ray energies and allow us to probe 
the high energy part of the HESS J1745-290 spectrum. The 
datasets are described in Table[T] Most of the data were taken in 



Table 1. Details of the observation of the GC region with 
H.E.S.S. for the 2004, 2005 and 2006 observation cam- 
paigns. 



Year 


range 


0/ 


-L obs^ 


N7I 






n 


n 


(hrs) 






2004 


0-60 


21.3 


28.5 


1516 


35.6 


2005 


0-70 


26.6 


51.7 


2062 


41.2 


2006 


0-50 


19.1 


12.7 


607 


23.0 


All 


0-70 


23.0 


92.9 


4185 


60.7 



" 9^ denotes the zenith angle 
' Live time of the data sample 
Number of detected 7-rays after background substraction 
Significance of the 7-ray excess N-, calculated according to 
the prescription of Li and Ma (|1983|l 



"wobble mode" where the telescope pointing is typically shifted 
by ± 0.7° from the nominal target position. The dataset used for 
the analysis was selected using the standard quality criteria, ex- 
cluding runs taken under bad and variable weather conditions, 
which would lead to large systematic errors. An additional qual- 
ity selection was applied using the cosmic ray rate during data 
collection. Runs with rates deviating by more than 3(t from the 
cosmic ray rate averaged over the three years of observation were 
removed in the subsequent analysis. After the run selection pro- 
cedure, the dataset amounts to 93 hours of live time. 
Two independent techniques are commonly used to select 7-ray 
events and derive energy spectra and lightcurves. The first tech- 
nique computes the "Hillas geometrical moments" of the shower 
image (Aharonian et al. I2005P and the second one is based on 
a semi-analytical model of showers, which predicts the expected 
intensity in each pixel of the camera (de Naurois et al. l2003|l . The 
shower direction, the impact point and the primary particle en- 
ergy are then derived with a likelihood fit of the model to match 
the data. Both analysis techniques provide an energy resolution 
of 15-20% with an angular resolution better than 0.1° above the 
analysis energy threshold. Results described in this paper are 
obtained with a combination of these two analysis techniques, 
the so-called combined Hillas/Model analysis (de Naurois et al. 
2003), which yields an improved background rejection, and were 
cross-checked against either technique alone. The background is 
calculated at each position in the field of view using the accep- 
tance corrected background rate from an annulus around that 
position (the OFF-source region) . More details on this so-called 
ring background method are given in Berge et al. (|2004| . The 
OFF-source regions do not overlap with known TeV sources in 
the GC field of view (Aharonian et al. I2006b|) and areas of dif- 
fuse emission (Aharonian et al. I2006all . 

The analysis of the whole 2004-2006 dataset shows an excess 
above the background of 4185 7-events in a 0.11°radius region 
centered on the GC (Table[l|. The total significance, calculated 
according to the method of Li and Ma p983p . is 60.7(t. Diffuse 
7-ray emission along the galactic plane is visible as shown by 
the distribution in 9, 9 being the angle of a reconstructed 7-ray 
relative to Sgr A* (Fig. [TJ excess above the background level 
outside the ON source region). A discussion in Aharonian et al. 
I|2006ap gives a possible interpretation of the emission as cosmic 
ray interactions in the central molecular zone (CMZ). A point- 
source -I- linear background fit to the 9 distribution has been 
performed to estimate the diffuse emission contribution to the 
signal. Extrapolating the linear component inside the ON-source 
region gives a diffuse emission contribution of 13% ± 1% within 
0.11°. Diffuse 7-ray emission located in the ON-source region 
was not subtracted in the following analysis, since its spectral 
feature (a simple power law with the same spectral index as the 
central source) would not influence the result (Aharonian et al. 
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l20U6ap . Any deviation from a power law spectrum or variability 
in the lightcurve would be solely due to the point source. 



(Eq. O and a smoothed broken power law (Eq.|3]): 
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Fig. 1. Distribution of the number of 7-rays derived with 
the combined Hillas/Model analysis as a function of 0, 
where is the angle between the 7-ray direction and the GC 
position. The dashed vertical line delimits the ON-source 
region. The horizontal line shows the expected number of 
cosmic ray background events. 



3. Energy spectrum of HESS J1745-290 

3.1. Spectral reconstruction 

The energy spectrum of HESS J1745— 290 has been derived fol- 
lowing the forward folding method. The forward folding method 
is based on a spectral shape hypothesis for the source spectrum 
(Piron et al. 120011 Djannati-Atai et al. 11999^ . Spectral points 
are then calculated using the adjusted spectral shape and the 
la error bars are computed using the error matrix of the fit pro- 
cedure. Relative systematic errors on the reconstructed spectral 
indices coming from the presence of broken pixels in the camera 
are less than 5%, whereas those coming from variations of the at- 
mospheric conditions are negligible. Thus, systematic errors on 
the spectral indices derived in the following analysis are taken to 
be 5%. Systematic errors on the integrated fluxes mainly come 
from the variations of the atmospheric conditions and the abso- 
lute calibration of the response of the telescopes. They amount 
to 10 to 20% (Aharonian et al. I2006dp . The systematic errors 
on the integrated fluxes above 1 TeV are taken to be 20%. A 
Monte-Carlo study of the reconstruction of the cut-off energy 
revealed a systematic bias linearly increasing with the cut-off 
energy: 



Ecut = (0.92 ± 0.01) X Ecut,true + (0.25 ± 0.05) TeV. (1) 



The systematic errors on the reconstruction of the cut-off energy 
amount to 17%. 



3.2. Results 

The data taken in 2004, 2005 and 2006 were compared to the 
folding of three distributions with the detector response: a power 
law (Eq. [5)), a power law with a high energy exponential cut-off 



dN 

Ie 

dN 

Ie 

dN 

Ie 



*°"(tIv)" 

E 



<l>o X 
<l>o X ( 



ITeV/ 

E \-ri 

VlTeV/ 



(2) 
(3) 
(4) 



where $o is the flux normalisation in TeV~ cm~ s^ , Ti the 
spectral indices. Ecut is the cut-off energy in Eq. [S] and Ebrcak 
is the break energy in Eq. |4] In Eq. |3l /3 is the strength of the 
cut-off. Except at the end of paragraph 3.2.1., I3 is taken equal 
to one. 

The measured spectrum for the whole three-year dataset ranges 
from 160 GeV, the energy threshold of the analysis, to 70 TeV 
(Fig. [2]). For the flrst time, with additional statistics, a devia- 
tion from a pure power law starts to be visible. The spectrum 
is well described by either Eq. |3] (equivalent of 23/26 d.o.f.) 
or Eq. m (equivalent of 20/19 d.o.f.). Fig.Elshows the HESS 
J1745-290 spectra with fits to Eq. [3] and Eq. g] The power 
law with an exponential cut-off fit yields "l>o = (2.55 ± 0.06stat 
± 0.40syst) X 10"^^ TeV"^ cm'^ s'S F = 2.10 ± 0.04stat ± 
O.lOsyst, a cut-off energy Ecut = (15.7 ± 3.4 

st at i 2 . 5syst ) TeV 

and an integrated fiux above 1 TeV of (1.99 ± 0.09stat ± 0.40syst) 
X 10~^^cm~^ s~\ The broken power law fit yields 4>o = (2.57 
± 0.07stat ± 0.40syst) X 10"^^ TeV'^ cm^^ s~\ Fi = 2.02 ± 
O.OSstat ± O.lOsyst, Fa = 2.63 ± 0.14stat ± O.lOeyst, a break en- 
ergy Ebreak = (2.57 ± 0.19stat ± 0.44syst) TeV and an integrated 
flux above 1 TeV of (1.98 ± 0.38stat ± 0.40syst) x 10-^^j,jjj-2 
s~^. The values of the cut-off energy Ecut and break energy 
Ebrcak are those corrected for the systematic bias mentioned 
in the previous section. By comparison, a power law spectrum 
gives a worse equivalent of 64/27 d.o.f. and yields a flux 
normalisation of (2.40 ± 0.05stat ± 0.40syst) x 10"^^ TeV"^ 
cm~^ s"'^, a spectral index of 2.29 ± 0.02stat ± O.lOsyst with 
an integrated flux above 1 TeV of (1.87 ± 0.05stat ± 0.40syst) 
X 10~^ cm"'^ s~^. The effect of introducing a high cut-off en- 
ergy in the power law spectrum does not change the integrated 
flux by more than 10%, less than systematic errors. Integrated 
fluxes found either with a power law with an exponential cut-off 
shape or a smoothed broken power law one are consistent with 
the values given in the previously published H.E.S.S. analysis 
(Aharonian et al. 120060(1 . The spectral parameters for the two 
curved spectral shapes, for year-wise data, are given in the next 
two subsections. 



3.2.1. Power law with an exponential cut-off 

Spectral parameters for the different years (differential flux nor- 
malisation $0, spectral index F, cut-off energy and integrated 
flux above 1 TeV, I(>1 TeV)) are summarized in Tabled The 
cut-off energy corrected for the systematic bias mentioned in 
section 3.1 is also shown. The energy ranges of the 2004 and 
2006 spectra do not allow the determination of a significant cut- 
off. The 2004 and 2006 datasets contain fewer runs collected at 
high zenith angles than the 2005 one. High zenith angle obser- 
vations provide an increased effective area at high energies and 
contribute to the high energy part of the spectrum. The 2005 
spectral index value is smaller than the values for 2004 and 2006, 
because of the correlation between the reconstructed spectral in- 
dex F and the cut-off energy induced by the fit procedure. Fig. 
|3] shows a 2-D plot of the fitted photon index F against the 
cut-off energy for each year's dataset. It can be seen that the 
spectral parameter values are compatible. Moreover, when fix- 
ing the cut-off energy to the uncorrected cut-off energy of 14.7 
TeV (see Table [S]), the spectral fits give F = 2.14 ± 0.07stat ± 
O.lOsyst, F = 2.04 ± 0.04stat ± O.lOsyst and F = 2.11 ± 0.09stat 
± O.lOsyst for the 2004, 2005 and 2006 datasets, respectively. 



4 



F. Aharonian et al.: Spectrum and variability of HESS J1745-290 




S- 2r 

■q. -I 

X I _ 

0) 

I ■ 

-1 : 
-2^ 



Energy (TeV) 



10 



Energy (TeV) 




Energy (TeV) 



Fig. 2. HESS J1745-290 spectra derived with the combined HiUas/Model analysis for the whole H.E.S.S. GC dataset 
covering the three years 2004, 2005 and 2006. The shaded areas are the la confidence intervals for the power law with 
an exponential cut-off fit (left) and the smoothed broken power law fit (right). The last points represent 95% confidence 
level upper limits on the flux. The fit residuals corresponding to the respective fits are shown on the lower panels. 



Table 2. Spectral parameters values with their respective statistical errors for a power law shape with an exponential 
cut-off, with (3—1. The fits are performed in the [160 GeV - 70 TeV] energy range. 



Year 


(lO-i^TevXm-^s-i) 


r: 


Ecut^ 

(TeV) 


J-Jcut ,tru<^ 

(TeV) 


I(>1 TeV)^ 
(lO-^^cm-^s-i) 


x7dof 


2004 
2005 
2006 


2.40 ± 0.10 
2.56 ± 0.09 
2.35 ± 0.16 


2.20 ± 0.07 
1.94 ± 0.07 
2.16 ± 0.11 


20.70 ± 11.80 

9.09 ± 2.13 
32.90 ± 39.50 


22.20 ± 11.80 

9.61 ± 2.13 
35.50 ± 39.50 


1.81 ± 0.14 
2.09 ± 0.16 
1.88 ± 0.22 


25/26 
33/25 
17/23 


All 


2.55 ± 0.06 


2.10 ± 0.04 


14.70 ± 3.41 


15.70 ± 3.41 


1.99 ± 0.09 


23/26 



" Flux normalisation 

' Reconstructed spectral index 
Reconstructed cut-off energy without corrections for the systematic bias 
Reconstructed cut-off energy corrected for the systematic bias 
Integrated 7-ray flux above 1 TeV 



Spectral indices are thus compatible with each other. The spec- 
tral shape given by Eq. [3] with /3 — 0.5 is motivated by some 
shock acceleration scenarios (see discussion in section 5.). A fit 
over the whole three years with /3 = 0.5 has been performed. 
It gives a reasonable xVdof of 30/23, $0 = (2.55 ± 0.07stat 
± 0.40syst) X 10"^^ TeV"^ cm'^ s'^, F = 1.91 ± O.OSstat ± 
0.09syst, a low value for the cut-off energy (corrected from the 
systematic bias) of Ecut = (4.0 ± 1.9stat ± 0.7syst) TeV and an 
integrated flux above 1 TeV of (1.98 ± 0.04stat ± 0.40syat) x 
lO-^^cm-2 s-^ 



3.2.2. Smoothed broken power law 



Spectral parameters (differential flux normalisation <l>o , spectral 
indices Fi and break energy and integrated flux above 1 
TeV, I(>1 TeV)) are summarized in Tabled The break energy 
corrected for the systematic bias mentioned in section 3.1 is also 
given. Parameters show a large dispersion (at most at the 3a 
level) but remain compatible with each other. The larger value 
of F2 for the 2005 data may reflect a possible steepening of the 
spectrum in the very high energy part. 
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Table 3. Spectral parameters values with their respective statistical errors for a smoothed broken power law shape. 
The fits are performed in the [160 GeV - 70 TeV] energy range. 



Year 


(lO-i^Tev'w-^s-i) 




J- 2 


Ebrcal^ 

(TeV) 


"FT 

-L-'brcak.truf^ 

(TeV) 


I(>1 TeVi^ 
(lO-i^cm-^s-i) 




2004 
2005 
2006 


2.39 ± 0.12 
2.64 ± 0.10 
2.28 ± 0.18 


2.18 ± 0.13 
1.73 ± 0.10 
2.27 ± 0.24 


2.51 ± 0.23 
3.07 ± 0.26 
2.19 ± 0.28 


2.40 ± 0.28 
3.35 ± 0.37 
2.04 ± 0.38 


2.33 ± 0.28 
3.37 ± 0.37 
1.94 ± 0.38 


1.79 ± 0.58 
2.11 ± 0.60 
1.85 ± 0.87 


25/19 
28/18 
14/18 


All 


2.57 ± 0.07 


2.02 ± 0.08 


2.63 ± 0.14 


2.61 ± 0.19 


2.57 ± 0.19 


1.98 ± 0.38 


21/19 



" Flux normalisation 

' Reconstructed spectral index 
Reconstructed break energy without corrections for the systematic bias 
Reconstructed break energy corrected for the systematic bias 
Integrated 7-ray flux above 1 TeV 



used to correct the measured spectrum: 



m" 0.2 



0.15 



0.1 



0.05 




F{E)^Fo{E) X exp(-r(£)) 



(5) 



where F(E) represents the observed spectrum after attenuation, 
Fo(E) the intrinsic spectrum and r(E) the optical depth of the 
7-rays toward the GC as a function of the energy. The optical 
depth r(E) depends on the number density of the ISRF photons 
and on the pair production cross section. 

Typically, 1% and 10% of 7-rays are absorbed at 10 TeV and 
50 TeV, respectively, with an increase of the attenuation fac- 
tor at higher energies. As expected, the spectral parameters do 
not change significantly after applying the correction. Thus, the 
curvature of the spectrum does not seem to be caused by 7-ray 
absorption and the cut-off (Eq. [3]) or the spectral break (Eq. |4]) 
is intrinsic to the source spectrum. 



I I I I \ I I I I I I I I I I I L 

1.8 2 2.2 2.4 2.6 



Fig. 3. Contours of constant likelihood of the power law 
with an exponential cut-off fit, as a function of the spectral 
index F and the inverse cut-off energy 1/Ecut for the 2004, 
2005 and 2006 datasets. The 68%, 95% and 99.9% contours 
are shown. The triangles denote the best fit positions. 



3.2.3. Spectral variability 

A refined study of the variations of the spectral index with time 
over the whole three years has been carried out using time inter- 
vals of roughly 5 hours, comprising ten consecutive runs. Each 
data subset has been fitted independently with a power law 
shape. The spectral index light curve has 25 points. A fit to a 
constant gives a of 29/24 d.o.f. confirming that the spectral 
index did not change significantly over the three years. 

3.3. Correction for absorption of very high energy ^-rays 

A recent calculation of the Galactic interstellar radiation field 
(Porter and Strong 2006) has shown that the infra-red radiation 
field near the GC is considerably enhanced compared to what 
was previously thought. Thus, some attenuation of very high 
energy 7-rays might occur at TeV energies. The attenuation 
coefficient accounting for very high energy 7-ray absorption by 
e'*'e~ pair production on the interstellar radiation field (ISRF) 
toward the GC was derived in Zhang et al. (|2006p and can be 



4. Search for time variability 

As mentioned in the introduction, fiaring of Sgr A* and QPOs 
were detected in various bands such as X-rays and IR. In this 
section, the "run by run" and "night by night" light curves are 
used to search for any variablility in the GC source activity. 
The run by run light curves are obtained by dividing the data 
into 28 minutes intervals and by computing the integrated fluxes 
above 1 TeV for each of these time intervals. The integrated flux 
is computed from the results of a power law fit with exponential 
cut-off whose normalization was varied in the likelihood mini- 
mization while the spectral index and cut-off energy were fixed 
to the values obtained for the 2004-2006 dataset. The fiux nor- 
malisation is adjusted for each 28 minutes time slice. The inte- 
grated fluxes above a fixed energy of 1 TeV are then calculated. 
The night by night light curves are derived in the same way as 
for the 28 minutes light curves except that the time intervals 
comprise runs that were recorded in the same night. 

4.1. Light curves and flare sensitivity 

The run by run integrated flux light curves of the GC covering 
the 2004, 2005 and 2006 observation seasons are displayed in 
Fig. |4l The fit of a constant to the data taken over the whole 
three years gives a of 233/216 d.o.f. and does not reveal any 
variability on time scales longer than 28 minutes. 
Because of the large error bars of the light curve points im- 
plied by the low statistics, the H.E.S.S. signal is only sensitive 
to relatively large amplitude flares. The flare sensitivity was es- 
timated by adding an artiflcial Gaussian with variable duration 
at, a time of maximum amplification to and an amplification at 
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Fig. 5. Maximum amplification factor A for a 3a flare de- 
tection as a function of the flare duration. A is obtained by 
a Landau fit (see text). The error bars are the correspond- 
ing la variations, depending on the time of the assumed 
flare. The maximum amplification factor decreases with in- 
creasing flare duration. 



maximum A to tlie H.E.S.S. light curve (LC): 



LC„od(t) = LC(t) X 1 + A X e 



(6) 



When varying the time of maximum amplification to along the 
LC, a distribution of the A values compatible with a flare detec- 
tion at a given significance is obtained. A fit to this distribution 
by a Landau function gives the most probable value of A and 
its corresponding la variations. Fig.[5]shows the maximum am- 
plification factor A compatible with no fiare detection at the Scr 
confidence level as a function of the fiare duration. Long flares 
involve a larger number of LC points in the constant flt and 
then increase the over the number of d.o.f. much more easily, 
so that A decreases with increasing fiare duration. Fig. [5] shows 
that an increase of the flux by at least a factor of two is neces- 
sary to detect flares of hour time scales. An increase of the 7-ray 
flux of a factor of 2 or greater was excluded at a confldence level 
of 99% during a Chandra flare night by the H.E.S.S. collabora- 
tion (Aharonian et al. I2008p . which is fully consistent with this 
result. 



4.2. Search for QPOs 

Four oscillation frequencies ranging from 100 s to 2250 s have 
been observed in the X-ray light curve of Sgr A* (Aschenbach 
et al. 12004 1) . These frequencies, if confirmed, are likely to cor- 
respond to gravitational cyclic modes associated with the ac- 
cretion disk of Sgr A*. The occurence of these frequencies was 
searched for in the data. First, the coherence time of the disk 
precession is assumed to be less than 28 minutes. A Rayleigh test 
(de Jager et al. I1989P is then performed on photon arrival time 
distributions for continuous observations of 28 minutes. Each 
28 minutes observation gives a Rayleigh power spectrum. The 
Rayleigh power averaged over 2004-2006 data is shown in Fig. |6] 
as a function of the frequency. Error bars are estimated by cal- 
culating the variance of the Rayleigh power distribution plotted 
at the corresponding frequency. The probed frequencies range 
from 1/28 min~^ to the inverse of the average time spacing be- 
tween two consecutive events of 1.2 min~^. No Rayleigh power 
is significantly higher than expected for noise at any probed fre- 
quencies. No significant peaks are seen at the 100 s, 219 s, 700 
s and 1150 s periods observed in X-rays. 

In a second analysis, the coherence time of oscillations is as- 




frequency (mm) 



Fig. 6. Rayleigh power plotted as a function of the fre- 
quency. Rayleigh power is normalized such that a pure noise 
spectrum results in unit power. The dotted line shows the 
fit to a constant of the Rayleigh spectrum. The arrows de- 
note the 100 s, 219 s, 700 s and 1150 s periods observed in 
X-rays. 



sumed to be of the order of a few hours. The Fourier power dis- 
tribution using Lomb-Scargle periodograms (Scargle 1982) for 
each night of the dataset is then constructed. Data are binned 
into 5 minutes intervals. The Fourier power averaged over the 
2004-2006 data is displayed in Fig. [7] as a function of the fre- 
quency. Frequencies tested range from 10~^ min~^ to 0.1 min^^. 
No significant oscillation frequencies are detected. 



5. Discussion and conclusions 

A strong signal has been detected from the GC region with the 
H.E.S.S. instrument. An energy spectrum has been measured 
with a difi'erential spectrum well described either by a power law 
with slope r=2.10 ± 0.04stat ± O.lOsyst with an exponential cut- 
off at 15.7 ± 3.4stat ± 2.5syst TeV, or a smoothed broken power 
law with photon indices ri=2.02 ± 0.08stat ± O.lOsyst, r2=2.63 
± 0.14atat ± O.lOsyst and a break energy at 2.57 ± 0.19stat ± 
0.44syst TeV. An integrated flux above 1 TeV of (1.99 ± 0.09stat 
± 0.40syat) X 10~^^cm~^ s~^ is derived using the power law 
spectrum with an exponential cut-off fit. No indication for vari- 
ability, flaring or QPOs has been found in the H.E.S.S. data, 
suggesting a non-variable emission of the GC region in the VHE 
regime. 

Different mechanisms have been suggested to explain the broad- 
band spectrum of the GC (Fig. |8]). Firstly, the stochastic accel- 
eration of electrons interacting with the turbulent magnetic held 
in the vicinity of Sgr A*, as discussed by Liu et al. (|2006p . has 
been advocated to explain the millimeter and sub-millimeter 
emission. This model would also reproduce the IR and X-ray 
flaring (Atoyan and Dermer |2004|). In addition, it assumes that 
charged particles are accreted onto the black hole, and predicts 
the escape of protons from the accretion disk and their accel- 
eration (Liu et al. 2006J. These protons produce tt" mesons by 
inelastic collisions with the interstellar medium in the central 
star cluster of the Galaxy. 

The cut-off energy found in the 7-ray spectrum could reflect a 
cut-off Ecut.p in the primary proton spectrum. In that case, one 
would expect a cut-off in the 7-ray spectral shape at Ecut — 
Ecut,p/30. The measured value given in section 3.2.1., with the 
parameter 13 equal to 0.5, would correspond in this scenario to a 
low cut-off energy in the primary proton spectrum of roughly 100 
TeV. It would correspond to a larger value of Ecut,p of ^ 400 TeV 
ff/3 = l. 

Energy-dependent diffusion models of protons to the outside 
of the central few parsecs of the Milky Way (Aharonian and 
Neronov l2005b|) are alternative plausible mechanisms to explain 
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Fig. 4. Run by run light curves of HESS J1745-290. From top to bottom: (1) Data covering the whole 2004-2006 time 
period. (2), (3) and (4) represent data covering the 2004, 2005 and 2006 time periods, respectively. The dotted lines show 
the fit to a constant integrated flux of the 2004-2006 lightcurve. 
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Fig. 7. Top panel: [10~^ min~^ - 0.1 min"^] Lomb-Scargle 
periodogram of the H.E.S.S. Sgr A* light curve aver- 
aged over the 2004-2006 nights of observation. Bottom 
panel: Fourier power distribution derived from the aver- 
aged Lomb-Scargle periodogram. No significant peak is vis- 
ible from the Lomb-Scargle periodogram and the of the 
exponential fit to the Fourier power distribution is 72/55 
d.o.f. 
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Fig. 8. A composite spectral energy distribution of the GC 
source. The radio data are from Zylka et al. (1995) and 
Falcke et al. (1998). The IR data are from Hornstein et al. 
(|2002p and Ghez et al (|2004p . Chandra data points come 
from Baganoff et al. ([Midla nd [2005]) . IGR J1745-2901 ones 
from Belanger et al. (^0051) . The HESS J1745-290 spectral 
points are shown by filled triangles. 
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U.K. Particle Physics and Astronomy Research Council (PPARC), 
the IPNP of the Charles University, the Polish Ministry of Science 
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of Namibia. We appreciate the excellent work of the technical sup- 
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Saclay, and in Namibia in the construction and operation of the equip- 
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the TeV emission observed with the H.E.S.S. instrument. They 
would lead to a spectral break as in the measured spectrum due 
to competition between injection and escape of protons outside 
the vicinity of the GC. A similar mechanism would explain the 
diffuse emission detected by H.E.S.S. along the galactic plane 
(Aharonian et al. I2006a|l . 

The recent discovery of the G359.95— 0.04 pulsar wind nebulae 
(PWN) located at 8" from Sgr A* also provides interesting mod- 
els as discussed by Wang et al. I|2006p and Hinton and Aharonian 
H2007P to explain the steepening of the measured spectrum of 
HESS J1745— 290. Inverse Compton emission due to a popula- 
tion of electrons whose energies extend up to 100 TeV might 
be responsible for at least a fraction of the TeV emission. The 
PWN model of Wang et al. and Hinton and Aharonian would 
imply a constant flux with time since the time scale for global 
PWN changes is typically much longer than a few years (more 
like centuries to millennia). 

The absence of variability in the TeV data suggests that the 
emission mechanisms and emission regions differ from those in- 
voked for the variable IR and X-ray emission. Moreover, the 
above-mentioned models can both accomodate a cut-off in the 
7-ray energy spectrum and predict the absence of variability in 
the TeV emission. They are thus viable scenarios to explain the 
strong TeV signal detected by H.E.S.S. in the GC region. 
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